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THE  U.S.  ARMY  (BRL'S)  KINETIC  ENERGY  PENETRATOR 
PROBLEM:  ESTIMATING  THE  PROBABILITY  OF 
RESPONSE  FOR  A  GIVEN  STIMULUS 


Thomas  A.  Mazzuchi 
Nozer  D.  Singpurwalla 


The  crew  compartment  of  an  army  vehicle  is  protected  by  an  armor 
plate.  It  is  desired  to  test  the  strength  of  this  armor  plate  in  order 
to  assess  its  appropriateness  for  use  in  the  vehicle. 

A  specimen  of  the  plate  is  taken  and  projectiles  arc  fired  at 
different  points  on  the  plate  at  different  striking  velocities.  If  a 
projectile  penetrates  the  armor  it  is  said  to  have  defeated  the  armor. 
Our  goal  is  to  determine  the  relationship  between  the  striking  velocity 
and  the  probability  of  penetration.  Due  to  the  expensive  nature  of  all 
items  involved,  this  goal  must  be  achieved  with  a  minimum  amount  of 
testing.  A  Bayesian  approach  for  solving  this  problem  is  presented 
here  and  illustrated  using  some  real  data. 
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THE  U.S.  ARMY  (BRL'S)  KINETIC  ENERGY  PENETRATOR 
PROBLEM:  ESTIMATING  THE  PROBABILITY  OF 
RESPONSE  FOR  A  GIVEN  STIMULUS 


Thomas  A.  Mazzuchi 
Nozer  D.  Singpurwalla 

1.  STATEMENT  OF  THE  PROBLEM 

The  following  statement  of  the  problem  is  based  on  our  several 
discussions  with  Dr.  Robert  L.  Launer  of  the  Army  Research  Office,  Re¬ 
search  Triangle  Park,  North  Carolina,  and  Dr.  J.  Richard  Moore  of  the 
Ballistic  Research  Laboratory  (BRL) ,  Aberdeen  Proving  Ground,  Maryland. 

The  crew  compartment  of  an  army  vehicle  is  protected  by  a  certain 
kind  of  material  which  we  will  refer  to  as  an  "armor  plate."  It  is 
desired  to  test  the  strength  of  this  armor  plate  so  that  we  may  be  able 
to  assess  its  appropriateness  for  use  on  the  vehicle. 

In  order  to  do  this,  a  10'  x  10’  specimen  of  the  armor  plate  is 
taken,  and  a  projectile  is  fired  from  a  gun  which  is  aimed  at  different 
points  on  the  plate.  In  Figure  1.1  below,  we  indicate  a  possible  firing 
pattern  according  to  which  the  gun  is  aimed. 

Typically,  the  distance  between  the  muzzle  of  the  gun  and  the 
target  is  about  200  meters,  and  rhe  velocity  of  the  projectile,  measured 
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Figure  1.1 — Illustration  of  a  firing  pattern  of  a  gun. 

between  two  conveniently  located  points  between  the  gun  and  the  target 
is  about  5000  feet  per  second. 

The  projectile  is  known  as  the  "penetrator,"  and  the  outcome  of 
each  firing  is  described  by  a  binary  variable  which  takes  the  value  1  if 
the  penetrator  defeats  the  target,  and  the  value  0  if  the  penetrator 
fails  to  defeat  the  target.  The  penetrator  induces  a  stress  on  the  armor 
the  stress  is  a  function  of  two  quantities,  the  "striking  velocity"  and 
the  "angle  of  fire."  The  striking  velocity,  also  known  as  the  "stimulus, 
is  the  velocity  with  which  the  penetrator  strikes  the  armor,  whereas 
the  angle  of  fire  8  (indicated  in  Figure  1.2  below)  is  the  amount  by 
which  the  armor  plate  is  tilted. 


Line  of  Fire 


Figure  1.2 — Illustration  of  the  angle  of  fire. 
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Both  the  armor  specimen  and  the  penetrator  are  very  expensive  and 
thus  the  testing  has  to  be  kept  to  a  bare  minimum.  One  strategy  that 
has  been  adopted  is  to  fix  the  angle  of  fire,  say  at  0*  ,  and  then  to 
fire  the  penetrator  at  different  striking  velocities.  After  each  firing, 
a  record  is  made  of  whether  the  penetrator  defeated  the  target  or  not. 

It  is  assumed  that  the  striking  velocity  can  be  measured  without  any 


error. 


2.  GOALS,  OBJECTIVES,  AND  SOME  COMMENTS 
ON  CURRENT  APPROACHES 

Given  that  our  goal  is  to  be  able  to  assess  the  appropriateness 
of  the  armor  plate  for  use  on  a  vehicle,  our  objective  should  be  to  es¬ 
timate  the  relationship  between  the  striking  velocity  (the  stimulus)  and 
the  probability  of  penetration  (a  response  of  1) .  This  is  illustrated 
in  Figure  2.1,  wherein  it  is  assumed  that  the  probability  of  penetration 
is  a  nondecreasing  function  of  the  stimulus. 

The  situation  described  above  is  identical  to  the  one  encountered 
in  "bioassay  experiments,"  and  "low  dose  radiation  experiments,"  in  which 
the  relationship  mentioned  before  is  known  as  the  quantal  response  curve. 
The  dose  level  of  a  drug  is  the  stimulus,  and  interest  generally  centers 
around  V  ,  the  stimulu"  at  which  the  probability  of  response  is  .5  . 

•  J 

Since  it  is  possible  to  subject  more  than  one  animal  to  a  particular 
dose  level,  the  number  of  tests  at  each  value  of  the  stimulus  can  be 
more  than  one.  Furthermore,  tests  are  often  conducted  at  several  dose 
levels,  and  thus  the  large  sample  theory  which  typically  justifies  in¬ 
ference  from  bioassay  experiments  is  adequately  substantiated. 


Probability  of 
Penetration 


Quantal  Response 
Curve 


Stimulus,  V 


Figure  2.1 — Probability  of  penetration  vs.  stimulus. 
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Despite  these  conspicuous  differences  between  bioassay  experi¬ 
mentation  and  the  problem  described  here,  the  methodology  and  techniques 
of  the  former  have  been  directly  adopted  for  use  in  the  latter.  In  so 
doing,  a  serious  compromise  has  been  made — the  estimation  of  V  , 
rather  than  the  entire  quantal  response  curve,  has  been  made  the  domi¬ 
nant  issue  of  the  kinetic  energy  penetration  problem.  Specifically,  the 
BRL's  commonly  used  "Langley  Method"  [Rothman,  Alexander,  and  Zimmerman 
(1965,  pp.  55-58)]  and  the  "Up  and  Down  Method"  [op.  cit . ,  pp.  101-103] 
focus  exclusive  attention  on  the  estimation  of  V  . 

•  3 

The  typical  approach  used  in  bioassay  for  estimating  V  is  to 
assume  that  the  probability  of  response  p  is  an  arbitrary  nondecreasing 
function  of  the  stimulus  V  ,  specified  via  the  relationship 

p  =  F ( (v-p) / o)  , 

where  F  is  a  distribution  function  determined  by  a  symmetrical  density 
function  with  location  parameter  y  and  scale  parameter  a  .  Often  F 
is  taken  to  be  the  normal  distribution  function 

x  _A  2 

F(x)  =  /  -  e  ~s  ds  , 

-00  /2tt 

*“X  ”1 

or  the  logistic  distribution  function  F(x)  =  (1  -  e  ) 

The  data  from  a  bioassay  experiment  consists  of  n^  ,  the  number 
of  subjects  receiving  stimulus  ,  i=l,...,K  ,  and  ,  j*l,...,n^  , 

where 

=  1  ,  if  the  j  subject  responds  under  stimulus  V^,  and 
=*  0  ,  otherwise. 

Civen  the  data  *  i=,»*-*»K  ,  j=l,...,n^  ,  the  param¬ 

eters  y  and  0  are  estimated  using  the  method  of  maximum  likelihood, 
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under  the  assumption  that  the  test  results  can  be  judged  independent. 

Once  p  and  0  are  estimated,  the  estimation  of  V  follows  from 

•  -/ 

the  fact  that  F  ,  the  tolerance  distribution,  has  been  specified.  Non- 
parametric  and  robust  estimators  of  V  ,  ,  such  as  the  Spearman-Karber 
estimator,  the  L-estimator,  the  M-estimator,  and  the  Tukey  Biweight  es¬ 
timator,  have  also  been  obtained,  all  under  the  assumption  that  the 
density  function  giving  F  is  symmetric.  These  estimators  have  been 
discussed  by  Miller  and  Halpern  (1979).  Furthermore,  it  has  been  empir¬ 
ically  shown  that  for  the  estimation  of  V  ^  it  does  not  matter  what 
specific  form  is  chosen  for  F  ;  many  of  the  commonly  used  nonparametric 
estimators  yield  identical  estimates  of  V  .  ,  as  long  as  symmetry  is 

•  J 

assumed . 

A  drawback  of  the  assumption  of  symmetry  is  that  the  estimate  of 
the  probability  of  response  when  the  stimulus  is  zero  is  nonzero. 

Whereas  this  may  not  be  too  disturbing  in  bioassay  with  its  emphasis  on 
V  ^  ,  in  the  problem  considered  here  and  the  low  dose  radiation  experi¬ 
mentation,  such  an  estimate  would  be  clearly  unacceptable.  A  zero  value 
of  the  stimulus  should  correspond  to  a  zero  value  for  the  probability  of 
response. 

In  view  of  the  above  difficulty,  the  paucity  of  data  at  each 
level  of  the  stimulus,  and  our  inability  to  specify  a  functional  form 
of  F  which  has  some  practical  merit,  we  are  motivated  to  advocate  a 
Bayesian  approach  for  the  solution  of  this  problem.  Our  approach  is 
described  in  Section  3. 
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3.  AN  OUTLINE  OF  A  BAYESIAN  APPROACH 


A  Bayesian  approach  to  the  bioassay  problem  was  first  proposed  by 
Kraft  and  Van  Eeden  in  1964,  and  was  more  fully  developed  by  Ramsey  in 
1972.  We  consider  here  the  theme  proposed  by  Ramsey;  extensions  of 
this  theme  are  considered  by  Shaked  and  Singpurwalla  (1982) . 

Let  0  H  V.  <  V.  <  ...  <  V  <  V  <  00  ,  be  M  distinct  levels 
of  the  stimulus  at  which  the  target  (armor  plate)  is  tested;  M  is 
chosen  in  advance.  The  outcome  of  a  test  at  is  described  by  a 

binary  (0,1)  variable  ,  where  X^  =  1  if  the  penetrator  with  a 
striking  velocity  defeats  the  target.  Let  p^  =  PfX^l}  ,  i=l,...,M 

and  without  loss  of  generality,  we  assume  that 

0  =»0  <  P1  <  P2  <  •••  <  »M  <  Vl  =  1  ;  <3-1) 


it  is  always  possible  to  choose  V,  and  Vw  which  satisfy  the  above 

1  M 

inequality . 

Given  X  =  (X^,...,X^)  ,  one  goal  is  to  estimate  the  unknown 
P± * s  ,  i=l,...,M  ,  subject  to  the  inequalities  (3.1).  Another  goal  is 
to  estimate  p_.  ,  for  some  j^i  ,  such  that  if  c  ,  the 

estimates  satisfy  p^  <  p^  <  >  i=l,...,M  ;  this  pertains  to  esti¬ 

mating  the  probability  of  response  at  a  stimulus  where  no  target  was 
tested.  Yet  a  third  goal  would  be  to  estimate  the  largest  stimulus, 
say  ,  for  which  £  a  ,  where  0  <  a  <  1  is  specified. 

Ramsey's  approach  for  achieving  the  above  goals  is  to  assign  a 
Dirichlet  as  a  prior  distribution  for  the  successive  differences 

P^»  P2~Pi»  •••»  P^-P^  ^  .  and  then  to  use  the  modal  value  of  the 

t 
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resulting  joint  posterior  distribution  as  a  Bayes  point  estimate  of 

(p^,...,p^)  .  The  modal  value  is  computed  with  the  inequalities  (3.1) 

being  satisfied.  The  modal  value  of  the  posterior  distribution,  if 

unique,  is  also  known  as  the  generalized  maximum  likelihood  estimator 

[see  DeGroot  (1970,  p.  236)],  and  is  used  as  a  Bayes  estimator  when  we 

do  not  wish  to  specify  a  particular  loss  function.  Having  estimated 

the  p.'s  ,  the  estimation  of  p.  and  V  is  undertaken  via  an  inter- 
i  j  a 

polation  procedure. 

Specifically,  if  a  >  0  ,  i=l,...,M  ,  and  3  >  0  are  cor,:  .-.s 

rM+l 

such  that  =  1  >  then  the  prior  density  function  ir  is  oi 

form 


TT  « 


M+l 

n 

i=l 


a. 

i 


(3.2) 


It  is  important  to  note  that  when  averaging  according  to  tt  in- 

M  M+l 

tegration  must  be  done  with  respect  to  H.  ,  dp.  /II.  (p.  -  p.  ,)  . 

i=l  *i  i=l  Ki  1 l-l 

Since  M  has  been  prechosen,  the  stopping  rule  is  clearly  de¬ 
lineated,  and  so  the  likelihood  for  the  response  probabilities  at  the 
observed  stresses  is 


M  X.  1-X. 

n  p/  (1  -  P.)  1  •  (3.3) 

i=l 

The  joint  density  function  of  the  posterior  distribution  of 
P^,...,PM  is  proportional  to  the  product  of  the  prior  density  function 
(3.2)  and  the  likelihood  function  (3.3).  Thus 
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8 


(3.4) 


f(pi’  “•  ,PM  I  Xl»  * 
M+l  X . 

“  n  p  1  (i  -  p  ) 
i=l 


1_xi  r(B) 

M+l 

n  r<3a.) 

•  i  -L 


M+l 

n 

i=l 


-  P 


i-1 


Ramsey  has  not  been  able  to  obtain  the  posterior  marginal  dis¬ 
tributions  of  p^  ,  i=l,...,M  ,  nor  has  he  commented  on  any  aspects  of 
these  distributions.  He  uses  a  nonlinear  programming  algorithm  to  ob¬ 
tain  (p  , ...,p  )  ,  the  modal  value  of  (3.4),  subject  to  the  constraint 

that  p  <  p0  <  . . .  <  p  ;  this  is  his  Bayes  estimator  of  (p, ,...,p„)  • 

1  =  Z  =  =  M  1  M 

In  contrast  to  this  Mazzuchi  (1982)  has  been  able  to  obtain  all  the 
moments  of  the  marginal  posterior  distribution  of  the  p^  ,  i=l,...,M  . 
This  work  of  Mazzuchi’ s  represents  an  extension  of  Ramsey's  results, 
and  is  one  that  takes  us  a  step  closer  to  a  fully  Bayesian  analysis. 

The  moments  can  be  used  to  approximate  the  marginal  posterior  distri¬ 
butions  of  the  p/s  using  the  techniques  given  in  Elderton  and  Johnson 
(1969) .  The  approximated  posterior  distributions  give  us  a  measure  of 
uncertainty  associated  with  our  using  the  first  moment  of  the  marginal 
posterior  distribution  of  p^  ,  i=l,...,M  ,  as  our  Bayes  estimate  of 
p^  .  The  first  moment  of  the  marginal  posterior  distribution  is  used  as 
a  Bayes  estimator  when  we  are  willing  to  assume  the  square  error  as  a 
loss  function.  The  formulae  for  the  moments  and  their  use  for  approxi¬ 
mating  the  marginal  posterior  distributions  are  given  in  Appendix  A. 

The  computational  effort  required  to  compute  the  moments  men¬ 
tioned  above  increases  with  M  .  Thus  there  is  a  trade-off  between  the 
convenience  of  using  an  optimization  algorithm  to  obtain  the  modal  value 
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of  (3. A),  versus  the  laborious  computational  effort  involved  in  obtain¬ 
ing  several  moments  of  each  of  the  M  posterior  marginal  distributions. 
The  optimization  algorithm  cited  above  is  based  on  the  "Sequential 
Unconstrained  Minimization  Technique"  (SUMT)  of  Fiacco  and  McCormick 
(1968).  A  computer  code  which  adopts  SUMT  for  the  problem  considered 
here  is  described  by  Mazzuchi  and  Soyer  (1982).  This  code  can  also  be 
used  for  the  computation  of  the  moments  of  the  marginal  posterior  dis¬ 
tributions  of  the  p^  ,  i=l,...,M  . 


3.1  Specification  of  the  Prior  Parameters 

In  order  to  implement  the  Bayesian  procedure,  we  need  to  specify 
the  prior  parameters  ,  i=l,...,M  ,  and  8  ,  given  in  (3.2).  In  order 
to  do  this,  we  observe  (see  Ramsey)  that  ik  =  p^  -  ,  i=l,...,M  , 

has  a  beta  distribution  on  the  unit  interval  (denoted  as 
uA  ~  Beta(3ai,  3(1  -  a_L) ;  0,1))  , 

r(S)  ^i  B(l-a  ) 

f(ui;  ^i’  =  r(Ba.)  r(3(l-a. ) )  ui  (1-ui)  »  0  ==  ui  -  ] 


with 


E(u^)  =  ,  and 


(3.5) 


Var(u^)  = 


ai(l  -  ai) 


(3.6) 


(3+1)  • 

If  P*  denotes  our  best  prior  guess  about  ,  consistent  with 

the  fact  that  the  P*’s  increase  in  i  ,  then  the  a.'s  can  be  ob- 

i  l 

tained  via  (3.5)  as 


and 


a  =  p* 
1  1 


a.  =  P*  -  P*  ,  i=2 , . . . ,M  , 

1  i  i-l  ’  ’  * 


Vfl  ~  1  ■  PM  • 

In  order  to  choose  the  parameter  0  ,  we  need  to  have  some  idea 
about  the  uncertainty  associated  with  our  choice  of  p*  .  This  in 
practice  can  be  done  in  one  of  the  following  two  ways: 

(i)  Suppose  that  in  addition  to  p*  ,  our  best  guess  about  the 
variance  of  p^  is  Va-r(p^)  .  Then,  substituting  =  p* 

in  (3.6),  we  have 

a.  (1  -  a  ) 

Var(ui)  ■  ',ar(pi)  ■  imT  - 


so  that 


p*(l  -  p*) 


-  1  ,  if  0  >  0  , 


g  =  )  Var(p1) 

\0  ,  otherwise. 

Note  that  0=0  corresponds  to  the  case  of  isotonic 
regression . 

(ii)  Often  in  practice  [cf.  McDonald  (1979)],  associated  with 
the  best  guess  value  p*  ,  a  user  is  able  to  specify  two 
numbers  a’£  >  0  and  b*  <  1  ,  such  that  for  some  y^ 
(specified  by  the  user),  0  <  y^  <  1  , 

P(a*  <  PL  c  b*)  =  1  -  yx  . 

Since  p^  -  Beta(0ou,  0(.l-ou) ;  0,1)  ,  given  p*  ,  we  set 
=  p*  ,  and  find  that  value  of  0  such  that 
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\  r(6a.)  r(0(l-a  ))  P1  (1'pl)  dpl  '1 

a*  1  1 

(3.7) 

Suppose,  further,  that  for  any  one  or  more  of  the  indices 
i  ,  i=2,...,M  ,  a  user  is  also  able  to  specify  two  numbers 
a*  >  pt  i  ,  and  b*  <  1  ,  such  that  for  some  y^  (speci¬ 
fied  by  the  user),  0  <  y^  <  1  , 

P(a*  <  (p*  |  p*_1,...,p*)  <  b*)  =  1  -  Y±  . 

Then,  using  the  fact  (see  Ramsey)  that 

(pi  I  P*^)  -  Beta(Bai#  6(1  -  ^  -  ...  -  c^);  p*^,  1) 

=  f(pi  I  P J_1 ;  B,“i)  »  say, 

we  can  find  the  smallest  value  of  6  ,  6*  ,  which  satisfies 
(3.7)  and  (3.8),  where 

bi 

f  f(pi  I  p£  6, a .)dpi  =  l  -  y.  ,  (3.8) 

ai 

with  a±  =  p*  -  p*_x  ,  i=2 . M  . 

A  computer  code  which  determines  the  smallest  value  of  6  described 
above  is  available;  the  details  of  this  program  are  given  by  Mazzuchi 
and  Soyer  (1982).  Our  reason  for  choosing  the  smallest  value  of  6 
stems  from  the  fact  that  large  values  of  6  give  a  very  strong  prior, 
with  the  result  that  even  a  large  amount  of  failure  data  will  not  change 


our  prior  distribution. 


3.2  Interpolation  Procedure  and  the 
Estimation  of  Quantiles 


Let  the  M-dimensional  point 


(P^»  •  •  •  >Pj.j) 


(p  ,...,p  )  ,  if  the  mode  of  the  joint  posterior  is 

^(p  ,...,p  )  ,  if  the  first  moments  of  the  marginal  poste- 
k  1  M 

nor  are 


used  as  the  Bayes  estimator  of  (p  , ...,pM)  . 

Suppose  that  we  wish  to  estimate  p^  ,  for  some  j^  i  ,  1=1,... ,M  , 
where  .  Let  p*  be  our  best  prior  guess  of  p  ,  the 

probability  of  response  at  a  nonexperimental  impulse  .  Then,  follow¬ 

ing  Ramsey ,  we  pick  p*  in  such  a  manner  that 


Pi+1  '  Pj 


Pi+1  Pj 


p*  -  p* 

1J _ li 

+  +  * 
p.  -  p. 

J  i 


(3.9) 


For  the  estimation  of  Va  ,  the  ath  quantile  (0  <  a  <  1)  ,  we 

first  see  if  there  is  an  observation  stimulus,  say  ,  for  which 

p*  =  a  .  If  so,  then  V.  is  our  Bayes  estimate  of  V  .  If  not,  we 
’  i  J  a  ’ 

determine  the  pair  of  observational  impulses,  say  and  »  f°r 

which  p^  <  a  <  P^+^  •  Since  the  probability  of  response  curve  is  as¬ 
sumed  to  be  increasing,  the  straight  line  segment  joining  the  points 
4  4  4 

0,  p^,  ...,  p^,  Pi+^>  . ..,  PM>  1  >  will  be  an  increasing  function  of  i  . 
We  shall  find  that  value  of  the  impulse,  say  V*  ,  <  V+  <  ,  for 

which  p*  =  a  . 


AO 


4.  APPLICATION  TO  SOME  BRL  DATA 


In  Appendix  B  we  present  eight  sets  of  data  labelled  1,  2,  3,  4, 
6,  7,  8,  and  9,  pertaining  to  60  kinetic  energy  penetration  tests. 

These  data  were  given  to  us  by  Dr.  Moore  of  BRL  and  have  been  carefully 
sanitized  to  maintain  confidentiality.  Data  sets  labelled  5  and  10, 
also  given  to  us  by  Dr.  Moore,  have  been  eliminated  from  consideration 
because  the  striking  velocity  for  these  data  is  much  too  different  from 
those  of  the  other  sets.  All  the  10  sets  of  data  were  obtained  sequen¬ 
tially  over  time,  in  the  sense  that  data  set  1  was  the  first  one  to  be 
obtained,  followed  by  data  set  2  (obtained  after  some  lapse  of  time), 
and  so  on,  until  we  reach  data  set  9,  which  is  the  last  considered  here. 
To  the  best  of  our  knowledge,  all  eight  data  sets  are  assumed  to  have 
been  collected  under  identical  conditions.  That  is,  there  is  no  indica¬ 
tion  that,  except  for  differences  in  striking  velocity,  the  material  and 
the  methods  of  testing  used  for  data  set  1  are  different  from  those  used 
in  data  set  2,  and  so  on.  This,  plus  the  sequential  nature  of  the  data, 
enables  us  to  use  the  posterior  obtained  from  one  data  set  as  the  prior 
for  the  next  set,  and  so  on,  until  we  obtain  the  posterior  using  data 
set  9,  which  then  gives  our  final  estimate  of  the  response  curve. 

Data  set  1  consists  of  13  observations  taken  at  striking  veloci¬ 
ties  ranging  from  128.60  (in  some  unspecified  units)  to  166.16.  The 
result  of  each  test  is  indicated  by  a  binary  variable  .  The  best 
prior  guess  values  p*  ,  necessary  to  choose  the  prior  parameters  ou  , 
were  not  specified  by  BRL.  However,  what  appears  to  be  reasonable  is  to 
assume  that  the  probability  of  response  at  a  striking  velocity  of  100  is 
close  to  zero,  and  that  at  a  striking  velocity  of  200  it  is  almost  1. 


Thus  we  make  an  arbitrary  choice  for  p*  ,  say  p*q  ,  by  letting 
P^q  =  1  -  exp[-.07(V^  -  100)3  .  Data  on  striking  velocities  outside  the 
range  of  100  to  200  were  excluded.  Despite  this  arbitrary  choice  of 
p^o  ,  we  shall  see  how  even  a  scant  amount  of  data  significantly  changes 
the  posterior  response  curve,  provided  that  the  smoothing  parameter  3 
is  not  too  large.  Three  values  of  3  were  also  chosen  arbitrarily; 
these  are  1,  10,  and  25.  Recall  that  small  values  of  3  tend  to  em¬ 
phasize  the  data,  whereas  large  values  of  3  tend  to  emphasize  the 
prior  distribution.  In  Appendix  B  we  show  our  analysis  for  the  case 
of  3  =  10  . 

Since,  in  reality,  the  data  are  generated  sequentially  over  time, 
our  first  step  would  be  to  revise  the  best  prior  guess  values  p*^  , 
i^lj.-.jbl  ,  based  on  data  set  1  alone.  The  posterior  (modal)  values 
corresponding  to  the  striking  velocities  of  data  set  1,  pT^  ,  will  be 
the  revised  values  of  p*^  ,  for  i=l,...,13  ;  these  are  given  in  column 
5  of  the  table  in  Appendix  B.  The  revised  values  of  p*^  ,  for  i=4,...,61  , 
are  obtained  via  the  interpolation  formula  (3.9),  using  pt^  ,  i=l,...,13  , 
and  p*Q  ,  1=14 ,...,61  .  Let  the  revised  values  of  p*^  ,  i=l4,...,61  , 
be  denoted  by  p*^  ;  these  too  are  shown  in  column  5  of  the  table  in 
Appendix  B. 

Upon  receiving  data  set  2,  we  revise  the  values  p*^  ,  i=l4,...,19  , 
by  the  posterior  modal  values  corresponding  to  the  six  striking  veloci¬ 
ties  of  data  set  2.  We  denote  these  revised  values  by  p^  *  i=14,...,19  ; 
these  are  given  in  column  5  of  the  table  in  Appendix  B.  The  revised 
values  of  p*^  ,  i=20,...,61  ,  arc  obtained  by  interpolation,  using 
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,  i=l, ... ,13  ,  p^2  ,  i=l4 , 19  ,  and  p*x  ,  i=20,...,61  ;  we 
denote  these  revised  values  by  p^  >  i=20,...,6l  ,  and  show  them  in 
column  6. 

We  continue  the  above  scheme  of  systematically  revising  the 
p^#'s  ,  either  via  the  posterior  modal  values  or  by  interpolation,  until 
we  incorporate  the  effect  of  all  eight  sets  of  data.  Data  set  9,  the 
last  one  considered  here,  consists  of  eight  observations  taken  at  start¬ 
ing  velocities  ranging  from  =  144.83  to  =  198.94  .  The  pos¬ 

terior  modal  values  corresponding  to  the  striking  velocities  of  data  set 
9,  Pi8  ,  i  =54,..., 61  ,  are  given  in  column  12;  the  interpolated  values 
p*2  required  to  obtain  the  p^g's  are  given  in  column  11.  Since  the 
p*2's  incorporate  the  results  of  the  previous  seven  sets  of  data,  we 
claim  that  the  final  posterior  modal  values  p^g  ,  i=54,...,6l  ,  are 
based  on  the  results  of  all  the  testing.  Had  we  ignored  the  sequential 
nature  of  the  data  and  computed  the  posterior  modal  values  by  using 
Bayes  Theorem  on  the  best  prior  guess  values  p*^  ,  i=l,...,61  ,  then 

the  posterior  modal  values  corresponding  to  V..  through  V,,  would  be 

14  51 

different  from  the  p*#  values,  i=l4,...,61  ,  given  in  the  table.  This 
difference  is  due  to  the  interpolation  scheme  that  is  used  to  constantly 
revise  the  best  prior  guess  values,  when  we  consider  the  data  sets 
sequentially. 

A  plot  of  p^g  versus  ,  i=54,...,61  ,  represents  our  final 
estimate  of  the  quantal  response  curve.  Estimates  of  the  probabilities 
of  response  at  striking  velocities  different  from  ,  i=54,...,6l  , 
can  be  obtained  using  the  interpolation  formula  (3.9).  When  we  use 
the  interpolation  formula  to  obtain  an  estimate  of  p.  ,  for  some 
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j=l,...,53  ,  we  need  to  specify  a  value  p*  ,  the  best  prior  guess  value 
of  Pj  .  Suppose  that  the  index  j  appears  in  data  set  k  ,  for  some 
k  <  9  ;  then  for  p*  we  will  use  p+^  .  In  so  doing,  we  will  have  in¬ 
corporated  the  effect  of  the  last  data  set,  data  set  9,  in  our  obtaining 

the  estimate  of  p.  ,  and  thus  achieve  a  certain  amount  of  smoothness. 

J 

Note  that  the  effect  of  the  data  sets  between  k  and  9  is  already  pres¬ 
ent  in  our  estimates  pt^  ,  i=54,...,61  ,  and  these  are  used  in  our  in¬ 
terpolation  scheme.  For  example,  suppose  that  we  wish  to  estimate  the 
probability  of  response  at  a  striding  velocity  of  158.52.  This  striking 
velocity  occurs  in  data  set  2,  and  lies  between  the  striking  velocities 
148.97  and  159.15  of  data  set  9.  The  index  j  corresponding  to  the 
value  158.82  is  17.  To  use  (3.9),  we  identify  P*+^  and  Pi+i  as 
being  .70499  and  .53014,  respectively,  p*  and  p+  as  .62881  and 

.42386  (see  data  set  9),  and  p*  as  .64436  (see  data  set  2),  and  com- 

+ 

pute  Pj  as  our  estimate  of  p^  . 

In  Figures  4.1,  4.2,  and  4.3,  we  show  plots  of  our  Bayes  estimate 
of  the  probability  of  response  at  the  eight  striking  velocities  of  data 
set  9,  for  8  =  1,  10,  and  25  ,  respectively.  Also  shown  are  the  90% 
probability  of  coverage  intervals  for  each  estimate.  These  intervals 
are  obtained  using  the  moments  of  the  posterior  distributions  of  p^  , 
i=54,...,61  ,  and  then  using  the  techniques  of  Elderton  and  Johnson 
(1969)  to  approximate  the  posterior  distributions — see  Appendix  A.  On 
each  of  these  figures  we  also  show  a  graph  of  our  best  guess  values 


p*0  ,  i=l,...,61  ;  these  enable  us  to  see  how  the  data  have  changed  our 


prior  estimates.  We  observe  that  the  90%  probability  of  coverage  inter¬ 
vals  tend  to  be  small  in  the  middle  of  the  range  of  the  striking  velocities. 


probability  of  response  for 


HDNI  OX  S3MVOOS  O Z  -  NOIXD3S  SSOMO  -  C*9  ON 


f a •••■• X* a***aaa*aBa*aaaaaa*«a*a#**M#iaa 
aaa** aaaaa ■aaaaaaa** aaaa* aaaaa aaaaaaaaaaaaaaaaaaii aaaaaaa aaaa* 
■•■•••••■a aaaaaaaaaaaaaaaaaaaa aaaaa aaaaaaaaaaaaaaa ■■■■ 


•  Saaa 


■  . . 


••■■•••Ml 


■inimii 

aaaaaaaaaai 


iaiiamaataaWMaai 


I  «■■*>'-«**"*****"*****  **********  *****  ■■""■■■■■■■■■"■■■■■'•■■aaaaaaaa  aaaaa  aaaaa  aaaaa  aaaa  aaaai 
■ ■  •••iV.Vt***  aaaaa  aaaaa  •aaaaaaaa#  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaa 

■  . . .  feaaaa  a aaaa aaaaaaaa ■■  aaaaaaaaJ 

■  aaaat‘>«aaaaaaaaaa*aa  aaaaaaiaaaaaaaaaaaai  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  hiimm 

■  ■■■  art>a  aaaa  aaaa  aaaa  a  aa  aaaa  aaaa  aaaaaaaaaaaaaaaaaaaa  a  a  aaaaaaaa  aaaaaaaaaaaaaaaaaaaa  MtiiHlPMMIIIMIIPIP 
liHai/iaMmaa*aaaMaMMauHH«MnauiuaitanauiiaaBaaaaMHiMaaMamaaa*iMMMMMMiHM«aim(iMK: 
■ •aa*«r«aaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaaaaaaa aaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaa aaaC 
■ aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaiaaaaaiaaaaaaaaJ 


■a aaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaaaaaaaaaaaaaa aiiaaaaaaaaiaaaaiaaa iaaaaaallPiiSVRinSiiSSSViVli 
I •aaaaaaaaaaaaaaa* aaa aaaaa aaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaaaa a*a*aaaaaaaaaaaaaaaaaai*aaaa**aaa*saaaa* aaaaaaa*^ 
la«  «a»»aaai aaaaaaaaaa  aaaaaaaaaa iaaaa aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa<a aaaaaaaa# aaaaaaftaaa aaaa aaaaaa aaaaaaaa^ 


. . . 

MMMMHMiMHM«a*piailMflliat«IIIWMNA 

iaaaCJaMaaaiaaaaaaaaaaaaaa 
|MWiMM»anwwiii 
iia*  aiaaa  M«aa»aaaaataaaaaai 

- - - - n - laaaaaaaaaaaaaaaaa# - nr-i-nifTiirinuniiiiiiriiT-TiiirTiiimi  iia  t  ui  fciaaa  laaa  . . . 

|  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaiaaaaiiiaaaaaaiiKaaaaaitMaaaMaMaaMMiaii^WiiiiaiiiaKaMaiaMaaaaiaaH 

■  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aiaaa  aaaaaaaaaa  iaaaa  ••■aaaaan  aaaaaaaaaa  . .  aaaaaaaaaa  HaaiaaaM 

I  «ai*vva*a(  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaataaaaaa  aaaaaaaaaa  aaaaaaaaaa  if  aaaaaaa|Aaaaa  •>iaa«aaaaiaaaaaaaaaaaaa^^| 

■  aaaa<taaaaa  aaaaaaaaaa  aaaaaaaaaa  aiaaa  aaaaaaaaaiaaaaaaiaaaiaaa*  aaaaaaaaaa  aaaaaaaaaa  BaaaaiaaaiMaiiMHiMaaiMdraaaantaaiiamumaaaMaaaM 

!»■■■«-  -«•■■■  aaaaaaaaaa  . .  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  »i»««»aa^Baaaiif«M»aan«a»»aaaaaaaa«| 

liaHwam  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  iaaaa  aaaa  aaaaaaaaaa#  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  nawaadjitmaan  aaaaaaaaaaaaaaaaaaaa  B 

1  —  - - - - - - a# - — - - *"■ — ti  mi - tt — nraniaiifiB 

ji  aaarraaaaaaaaaaaaa*  aaaa  a  aaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaa *****  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaft»aaaa**aaa  aaaaaaaaaaaaaaaaaaaa  B 
l**!!^af!a<aaa*aisu*iaia*a*Maa«»iiaiMMMaaaHaaMaiMM«aa«aaia«t«uHana«»ntaMa«aaai«nu»n«Hii(iau««w'MM#«a>uu«ttMaaBMNl 

■  aaaa  —  >aaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa*  ■•aaaaMMaamarnmpa^iai  aaaaiaaaaaaaaaaaaaaa  I 

■  aaaaaaaaaa  aaaaaaaaaa  iaaaa  ■«»««  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  »aaMaaaM«aawaaa«aiMaaaaHi»aanaaatni»»aa  iXaaaWaa  aaaaaaaaaaaaaaa  aaa^J 

■  •aa*v*.aaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  a  a  aa#  ■  a  aaa  aaaa  •  aaaaaaaaaa  aaaaaaaa  ■  ■  aaaaaaaaaa  aaaa#  aaaaaaaaaa  aaaaa  aa»v?*Baat«wiaaa  a  aa  aaaa  ■  aaaaaaaaaa^! 

■  aaaar-iaaa  a  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaasaaaaaaaaal  aaaaaaaaaa  aaaaaaaaffi*aaaa~*aa(aaaaaaaaaaaaaaaaaa^H 

aa\-waaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaa  aaaaaaaaaa  nuaaaaataaMitaaaaaM^I 

■  aaa*  “  iaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaa  iaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaa  77vaaar*'*aai  aaaaaaaaaa  aaaaaaaaH 

■  •  ■»" -aaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaiaaa  iaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  iaaaaaa*XaMik(iaai  Ma«ita«*a»iBMaMH 
■•ataiawaa  aaaaaaaaaa  ■B«iaaa*ii  aaaaaaaaaa  aaaaaaaaaaaaiaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaMtofeaaC»<aaaaaa»s  aaaaaaaaaa  bmH 

- ■ - n - ruar - t  —  an«—  i»i  Timniair"i~nnnaaMH 

■  •  a>  i  a* -aa  aaa  a  aaaa  aaaaa  iaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaa  a  Maaa  aaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaa  aitia  mb  nnaaa  aaBMaS^I 

■  aaoaiaaaaa aaaaaaaaaaaaaaa aaaaaaaaaaaaaaaHaaaaaaaaiaaaaaauaaaaaaaam aaaaaaaaaa  aaa asaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaa  aaaaaaaaaaaaaaa  aaa^H 

■  •****•  aaaaaaaaaa  aaaa  aaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaaiaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaaaaaaa  aaa  as  «aanatMiinaaMaH 

■  aanai  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaraaaa  •  aaaa  aaaaa  aaaaa  aaaaa  aaaaa  iaaaa  a  aaaa*aaaaaaaaaaaa«#*aaa**aaaaasBaa**aB*aa*aa*a*Ba#aa*iaaaaa*a*a**a*^M 

■  ■  a<  in  iaaaa  aaaaaaaaaa  aaaaaaaaaa  aBaa^aaaa  aiaaa  aaaaaaaaaaaaaaa  aaaiiaaaaaaa  aaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaa  Him''  aaaaa  MMiaaaH 

■  ••>«ar«Baaa  ••aaa**aa*  aaa***«aaa  aaaavZaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaaaaaaaaaaat  »*ai»M«M»»f?faaaa««aaaaM»M 
|  ■  m>  iHaaaaa  aaaaaaaaaa  aaaaaaaaaa  a  aaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaai  aaaaaaaa  aaaaaaXaaaaaaaaaaaaa^H 

■  ■•  hi  aaaaaaaaaaaaaaa  ■aaaaaaaaaaaaar  ‘  aaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaa  »«aaaaaaaaaiaaM 

■  aaaaa  aaaaa  aaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaaiyaaaaaaaaaaaaaaaaataaaaaaaaaaa—  aaaaaaaaaaaaaaa  aaaaaaaa  aaaa— »a— «aaa— — M 

■  »a,iaa>a«aaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaa— aaaaaaaaaa  aX-aav»a*aaaaaaaaaaaiaaaaaaaa#aaaaaaaaaaaaaaaaaan  aaaaa  aaa  aaaa«a*aaa aaaaa aaaaa^M 

■  ••••■  aaa  aaaa  a  aaaaa  aaaaaaaaaa  aaaa*  a  aaaa  aaaaaaaaaaaaaaaaaaaa  •r.aMnaamaaaNBi  aaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaaa  aaaaa  aaaaa  aaaa;  aaaaMaaaaaaa^l 
|  ••••••aaaa  aaaaaaaaaa  aaaaa  iaaaa  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  a a^*M»aaaa aaaaaaaaaaaaaaaaaaaa  aaaaa  aaaaaaaaaa  MaMaaaaa  ati  aa  »aauaaaaa  aaiiaMaH 

■  ■  i'*B*ja  aaaa  aaaaaaaaaa  aaaaa  aaaaa  aaaa««  aaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaa  t.«aa*aaiaaaa*aa  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaa  aaa  aaaaaioaaaaaaaaaaaaia^H 

■  a*  la^aaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaa  a  aaaa  aaaaa  a  a  aaa  iMal  aaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaai  aaaaataaaaaaaM 
■•ataman#  aaaw  aaaaa  aaaaaaa— a  aaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaaa  >g»w>aaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaa  aaaaaaaaaa  aaaaa  aaaa-  a  aaaaaaaaaa aa^M 
B  —  iai.»am  aaaaa  aaaaa  aaaaaaaa— aaaar.  1 1<  at  aaaaa  aaaaa  aaaaa  aaaaa  fmwa— awaawa>«na«a»aa»  aaaaaaaaaa  aaaaa  aaaaaaaaaaaaaaa  aaa«»aaaaa  aaaaa  aaa^M 
■aamaaMaama— iai«aaaaaa«>aaa»aa^aa>ltaMaaaa»aaaat»Ma— awaiM—  aaaaa  aaaaa  i  aaaaaaaaaaaaaaa  aaa  a  aaaaaaaaaa  aaaaa  aaa  aaaaaaa  aaaaa  a— aaaai^M 

»»»aaaaaaa  aaaaaaaaaa  aaaaa  a  aaaa  aaaar -— •»  aaaaa  aaaaaaaaaaaaaaa  *-~r  ja"waa  aaaaaaaa#  aaa  aaaaaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaa  aaaaaaaaaaaaaaa  aaa  aaiaiJ 

■  •  a  ■  at*aa  awa  aaaaa  aaaaa  a  a  iaaaa  aaa  aaaaaaaaaa  naaaaaMiMaMMH||||MHM||M||MMMMMM|||M|||MMMMMHMMmmMMMMMHB|^H 


_  a  — .  aaaa aaa a.aa 
I  >>•«■»•■>  aaa  aaa 


•■a  aa •aaaaaaaaaa 

•  aaaaaaaa: - 

aa  aaaaa* 


U»aaaa.  aaaaaaaaaaaaaaaaaaaa  a  M*dPBa-H  aaaaaaa  aaaaa  aaaa  aaa  a»a  aaaaaaaa*  >  aaa  a  a  i  a  aaa  oJ 
- ar- * aaaaaaaaaaaaaaaaaaaa a«*aC= laaa aaaaaaaaaa aaaa* aaaa# aaaaaaaaaa aaaaaaaaaa H 

—^—tr**aaa*a*a  aaaaa  aaaaa  aaaaaaa  ■  aa  aaaaa  aaa  aa  aaaaa  aaaaa  aaaa#  aaaaa  aaaaa  aaaaa  aaaai 
I  mm  mm'  iaaaa  ’■<•••  mmmm  mw  >•  aa  ■■aaai  la-caai  aaaaa  aaaaa  aaaaa  aaaaa  a*Vaa  aaaaaaaaaa  aaaa#  aaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaaa  aaaaa  aaaaaaa! 

lniKiaaixaHM  «.  ««  •M,.ar'irl|uMi*iaa  aaaaa  aaaa*  a*7a*,vi*a*  aaaaa  aaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaa  aaaaa  aafafc*  aaaa*  alMIIIIRIIIINH 

1  aaaaaaaaaaaaaaaaaaaa  aaaaaaaaaa  a»  iaaaa  aa>»»a  *aa— wiaa»ia»a«M.  .aaaaaaaa  aaaa  aaaaa  ******************************  *afcii*a*g****^*B**a*****ia*BiB 

■  aaaaaaaaaa  aaaaa  ao  aaa  aaaaa  aaaaaaaaai- >**-.  »  ac  .an  aaaaa  aaaaa  aaaaaaaa  a*  »*»  *************  aaaaaaaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaa*  aaa  **B*ar?a***a  aaaaa  aaa  aal 

I  ■■*"  aaaaaaaaaaaaaaa  aaaai  aaaaa  m  Wanra»a»»»aaaa»aaa»Hti  ’  laaaaaa***************************  aaaaa  aaaaaaaaaa  aaa  aa  aaai#aaaaaa*a*aaaa  | 

■  aaaaa  aaaaaaaaaa  aaaaa  aaaaa  aaa—  *******r,^»*CL«a  aaaaaaaa  aa  aaaaa  ai^aamai  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aaaaa  aaaaa  aaaaa  aaa  aa  aaa  mm  *——■■■■» —J 
*■*•«  a  aaaa  •••••  aaaaa  a  *•■•■••■•  aaa  af , ,  *<r,Tl  •^V’aaaaaaaaaaaa  aaaaa  arVaa  AM  a  aaaaaaa  aaa  aaaaaaaaaa  •■■*•  aaaaa  aaaaaaaaaa  aaaaa  aaaaa  **»•-  I 
■  a  aaa  aaaaa  aaaaaaaaaa  aaaaaaaaaa  aaaaa  >  mm*  aa*  a*m  aa  aaaaaaaaaa  at— ^a  m  aaaa  aaaaaaaaaa  aaaaaaaaaa  a  aaaa  aaaaa  ■aaanaawaaaaa  aaaaa  aaayi 


m* aaaaaaaaai aaaaa aaaai aaaaa « 


aaaaaaaaaa^aaaaSaSal 
aaaaa* aaaaaaaaaa **i 


aaaaaaaa aaaaaaaaaa aaaaaaaa 


iBBaaBaaaaa«uaaaaaafe< 
aa aaa aa aaaa* aaaaa a 


RfRIR?!l!!!SVni^imniS9IR!9i*N««a«wVaaaiiBuaaBiM*aai 

. . .  '  -  -U3.  - /aaaaaaaaaa  aaaaaaarW  aaaaaaaaaa  aaaaa  aaaai— IP— ———I 

■HHBBaaaa  Iaaaa  aaaaaaa*  aaaaaaa  aaatoaaaaa-iaaaaaaaaaaiaaaaaaaa*  aaaaaaaaaaaaaaaaaaaa  laaaaaaaiaaaaaaaaaaaiaaaiaaaaia 
1 aaaaaaaaaa aaaaa aaaaa aaaaaaaaaa aaaaaaa aaaaaaaa aaaaaaaa aaa* /aaaaaaaaaa aaaaa ******************** aaaaa aaaaaaaaaa m 
|*a«ari»a  aaaaaaaa  aa,  aaa**  aaaaaaaaa'/aaaar-^***— •*****»*■  aaaa  «  aaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaaaaaaa’aaaaaaaaaaaaaaaa 
i  ••**>«#***  aaaaaaaaaa  aaaaaaaa*  aaaaa  ■■  ■■aaaa  ************************************************  a 

i  aaaaaaaaaa  aaaaa  aaaaa  aaaaa— aaaaa  aaa  aaaa*  *******>.!>***>■.**  a*  ***********************  ******************************  a 
■  aaaaaaaaaa  aaaaa  aaaaa  aaaaaaaaaa  aaaaa  aaaaaaa— aaaaa  a^^a—Kaia  aaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa—  a  aaaaaaa  aaaaaaaaaa  a 
■ aaaaaaaaaa aaaaa aaa aaaaaaa aaaaa aaaaaaaaaa aaaaaaaaaa aaaaaaaa aaaaaaaaaaaaaa aaaaaaaaaaaaaaa aaaaaaaa* aaaaaaaaaaaaaaa 

|a*aaaa**»a«a*ayaaa*aaaaaaaaaaaa*aaaaaaa*aaaaaaaa—— ———■■MBBB 


i#  aaaaaaaaaaaaaaaaaaaa  aaaaa  ■ 
a  aaaaa  ••aa~' aaaaa  aaaaaaa***  B 
4  aaaaa  aaa  ^aaaaa aaaaa aaa**f 


•  aaaaaaa  Bj 
a aaaaaaaa 
• aaaaa a*« 
iiaaaaa*# 
W  aaaaaaaa 


k  a  aaaaaaaaai 


- -  —  aaai - 

******  ’amaaat  aaaaa  aaaa*  aaaaa  aaaaaaa  - -  -  - 

—aaaaaaaaaaaaa  aaaaaaaaaaaaaaaaaaaa  aa  waimi885a>> 


idinMian  iikiMaai 
—  aaaa  aaaaa 


■  a  a 

■  fcaaaaaa* aaaaaaa  aaaaa  immH - , 

I— Mfw«ia«f  aaaaaaaaaaaaaaaaafTWWlMU  ■—■■■■  9S— 

yH**iii(ttma*(H*NM«M  ■••••••■»  a  aaaa*  a -a*a  at 

L*a«  ■m«an«iMt<  a  aaaaa  iimini  »  ••■  ■■  aaaaa  ai 


aaa# aaaaaaaaaa aaaaa 
tafti  aaaa  aaaaaaaaaaaaaaa 
kaaaaaiaai  aaaa  aaaaa 


aaaa  > aaaa a aaa aaai 
aaaSaaaSaBaaaaaai 


■■■aa  aaaaa  aaaaa  j 
aaaaaaaaaa*  aaaaa ■ 
aaaaaaaaaaaaaaaaa 


■  n»« 


■aaaaaaaaaaaaaaa 


aaaaa  aaaaaaaaaa  aaaaa  aaaai 
aaaaaaa aaaaa aaaaaaaaaaaaaaaaaaaa aaaaaaaa* 
>a*  ia  a  a  a  aaaai  aaaaa  aaaaa|MaaaMiaMaaaaaaa| 


■■■••■•a 

aaaaaaafl 

I  aaaaa  ad 
a  aaaaa aj 


■amiMaafll 

Bgssasssssl 


s«s  u 

m 


jaatni^P 


'Mim 

mm  in::! 


\ms 


nni 


finut 


Stuns 

aaaaa*** 


a*  aaaaa  awOTB 

paginal 

IBSSM 


%m\ 


46 


MDNt  OX  S3MvnOS  O Z  •  NOIXD3S  SSOdD  •  Cr 9  ON  ^T3avXlD 


In  Figure  4.4,  we  superimpose  the  plots  of  Figures  4.1,  4.2,  and  4.3,  in 
order  to  give  a  perspective  of  the  effect  of  3  in  our  computations. 

It  appears  that  our  Bayes  estimates  for  the  three  cases  of  3=1,  10, 
and  25  tend  to  converge  toward  each  other;  this  is  to  be  expected, 
since  we  have  61  observations  with  which  we  revise  our  prior 
probabilities . 
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APPENDIX  A 


/■ 


Moments  of  the  Marginal  Posterior  Distributions 


The  moments  of  the  posterior  distribution  of  ,  i=l,...,M  , 
have  been  obtained  by  Mazzuchi  (1982);  a  formula  for  obtaining  these  is 
given  below.  A  computer  code  which  facilitates  the  computation  of  the 
moments  is  described  by  Mazzuchi  and  Soyer  (1982) . 

Let  X  =  1  -  X  ,  i=l,...,M  ,  B(a,b)  =  T(a)r(b)  /  T(a+b)  ,  and 


•i  *M  r  ,  r.  M 

K  =  l  ...  I  (-1)  1=1  1  - 


rM 


rl=° 


r  =0 
M 


M  i  1 

n  B  l  X.  +  6a,  +  r.,  Ba  )  . 

t-i  l  A  J  33  1+1 1 


Then,  for  Jl=l,2,...  , 

,  i  X1  ^M 

E(pf)  =|  l  ...  I  (-1) 


rM 

ri 


r!=0 


rM=° 


M  I  i 

n  b  l 

i=l  \j  = 


^  X*  +  Ba.  +  r.,  Ba.+1 


) 


where 


X*  = 
J 


'X.  +  Z  , 
3 


'X. 

3 


3  =  s 
otherwise . 


These  moments  can  be  used  to  approximate  the  posterior  distribu¬ 
tion  of  p^  ,  f(p.j,)  ,  i=l,...,M  .  In  order  to  do  this,  we  consider  a 
system  of  frequency  curves  described  by  Elderton  and  Johnson  (1969) 
which  are  based  on  the  transforms  of  a  standard  normal  variate  Z  . 

The  system  of  curves  which  is  appropriate  to  our  problem  is  that  referred 
to  as  the  "bounded  system  of  curves,"  denoted  by  Elderton  and  Johnson 
(1969,  p.  123)  as  S  ,  and  described  by 

D 


Z 


Y  +  <5  Ln((p.  -  e)  /  (e  +  X  -  p.)]  ,  e  <  p.  <  e+X  , 


where  y  ,  6  ,  X  ,  and  £  are  parameters  whose  values  are  determined 
by  the  first  four  moments  of  f(p^)  about  its  mean. 

Hill,  Hill,  and  Holder  (1976)  give  a  computer  code  which  deter¬ 
mines  y  ,  6  ,  X  ,  and  £  from  the  first  four  moments  of  f(p^)  about 
its  mean.  Since  it  was  assumed  that  p.  ,  <  p.  <  p.,,  ,  we  estimate 
X  and  e  from  the  Bayesian  estimates  of  the  p^  ;  y  and  6  are  ob¬ 
tained  from  the  computer  code.  Having  obtained  these  parameters,  the 
distribution  f(p^)  is  obtained  from  Elderton  and  Johnson  (1969,  p. 
130)  as 


f(p±) 


V)]  ex!>[-  i  (T  +  ]  ’ 


£  <  p.  <  £  +  X  , 
*1 


where  N  in  our  case  is  1. 

In  order  to  obtain  the  approximate  (l-y)%  probability  of  coverage 
intervals  for  each  p^  ,  which  contain  its  Bayes  estimate  p^  (mode  or 
mean) ,  we  use  the  fact  that  since 

z  =  y  +  <S£n[(pi-£)/(£+X-pi)]  ,  £  <  p£  <  £  +  X  , 

Pt  “  X  exp[(^2.)  +  l]  +  £  . 

Thus,  to  find  two  numbers,  a  and  b  ,  such  that 

P{p*-a  ^  pi  £  P*+b}  =  1  -  5  , 


we  use 


P  -6Zn 


Y  <  2  <  n 


l~r~-  A  + 

\pi+b-£  I 


=  1  -  6  , 


and  solve  for  a  and  b  by  setting 


where  z.  ,,  is  the  (l-(6/2))th  percentile  of  a  standard  normal  distri- 
6/2 

bution.  Taking  c  =  max(a,b)  ,  we  form  our  interval 

Pr{p*  -  c<Pi<p^  +  c)  >_  1-6. 

These  intervals  may  not  be  symmetric  about  the  mean  or  modal  estimate. 
This  case  arises  when  the  boundaries  of  the  probability  of  coverage 
interval  exceed  the  boundary  of  the  variable.  In  such  cases  the  variable 
boundary  is  used  as  the  boundary  of  the  probability  of  coverage  interval. 
The  probability  of  any  symmetric  interval  about  the  mean  or  modal  esti¬ 
mate  may  be  obtained  by  proceeding  in  the  reverse  or  the  above  and 
evaluating  the  interval  for  the  standard  normal  variate. 


APPENDIX  B 


In  the  table  below  we  give  values  of  the  striking  velocity  V  , 
the  response  ,  and  the  best  prior  guess  values  p*g  ,  i=l,.-.,61  , 
for  the  eight  sets  of  data  described  in  Section  4.  We  also  show,  for 
8  ■  10  ,  the  revised  values  of  p|q  ,  p^  ,  or  p*^  based  on  data  set 
j  ,  j-1,2,3,4,6,7,8,9  . 
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